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A finite element simulation of sound attenuation in a 

FINITE DUCT WITH A PERIPHERALLY VARIABLE LINER 

by 

Willie R. Watson 

NASA-Langley Research Center 
Hampton, VA 

Using multimodal analysis, a variational finite element method is 
presented for analyzing sound attenuation in a three-dimensional finite 
duct with a peripherally variable liner in the absence of flow. A rect- 
angular element, with cubic shape functions, is employed. Excellent 
comparison between exact results and the finite element method is obtained 
in cases where an analytical solution exists. This study indicates that, 
once a small portion of a peripheral liner is removed, the attenuation 
rate near the frequency where maximum attenuation occurs drops significantly. 
Also, it was observed that the positioning of the liner segments affects 
the attenuation characteristics of the liner and that effects of the duct 
termination are important in the low frequency ranges. In general, the 
results indicate that the main effect of peripheral variation of the liner 
is a broadening of the attenuation characteristics in the mid-frequency 
range. Finally, it is concluded that, due to matrix size limitations of 
the presently available computer program, the eigenvalue equations should 
be solved out of core in order to handle realistic sources. 
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INTRODUCTION 

Lining the interior surfaces of aircraft engine ducts with acoustic 
treatment is a well established method for reducing internally generated 
aircraft engine noise. Large conimerlcal aircraft such as the L-1011, DC-10 
and Boeing 747 have successfully passed FAA certifications due, in part, to 
this concept and are positive evidence that acoustic liners can reduce air- 
craft noise effectively. However, increasingly stringent noise reduction 
goals require that these acoustic suppression techniques be continually 
refined and updated. Aircraft companies are now more than ever before in 
need of new ideas and methods to enhance the present state of liner 
technology. 

Initially, liner research was centered around uniform liners (ref. 1). 
Later, Zorumski and Lansing (ref. 2) realized that liners could be made more 
effective by taking advantage of impedance changes in axial segments. In 
light of this development, Zorumski (ref. 3) developed a theory to compute 
the attenuation in axially segmented circular and annular ducts. Several 
other investigators (ref. 4, 5, 6, 7 and 8) have since investigated axially 
segmented duct liners and their practical application. 

The present analysis was motivated, primarily, by the success of axially 
segmented liners and the desire for further effective methods for minimizing 
internally generated aircraft noise. In this work, a new type of liner 
variation, the peripherally variable liner, is investigated. In this type 
of liner, the impedance of the liner is allowed to vary around the duct 
perimeter, but remains constant in the axial direction. The impedance 
boundary conditions for a peripherally variable liner become boundary 
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conditions with variable coefficients and an exact analytical solution for 
the acoustic field cannot be determined. 

In this work, a finite element method is employed to extract the acoustic 
field and calculate sound attenuation in a three-dimensional rectangular duct 
with a peripherally variable liner. First, the governing Helmholtz equation 
and impedance boundary condition for zero mean flow are transformed into a 
single function, II, which has the governing Helmholtz equation and impedance 
boundary conditions as its stationary conditions. Next, the f-’*’'-i 2 element 
method is applied to approximate the functional II. Cubic functions which 
insure continuity of acoustic pressure and particle velocity throughout the 
duct are employed as shape functions. The requirement that the functional 
II be stationary results in a set of matrix equations which are solved to 
obtain the acoustic modes. Source and termination effects due to finite 
duct length are also included. 


ANALYSIS 

In this section, the mathematical expressions necessary for evaluation 
of the attenuation of peripherally variable liners will be developed. 

Statement of the Problem 

The duct to be analyzed and the Cartesian coordinate system to be used 
are shown in figure 1. The length, height, and width of the duct are d, h, 
and i respectively. The duct walls are acoustically lined with a specific 
acoustic admittance, 3(x,y}, which is independent of the axial coordinate 
z, but is a function of position along the duct perimeter so that: 


( 1 ) 


’3(x,oT 


Bi(xf 

3(x,Jl) 


02(x) 

S(0,y) 


03(y) 

B(h,yl 


R..(y)^ 


A constant exit admittance, 3g. is specified at the exit plane (z = d), and 
an acoustic pressure distribution, G{x,y), is input at the entrance plane, 

(z = 0). 

In separable geometries with constant values of 0 on the duct perimeter 
the solution to the acoustic field is given in terms of known trigometric 
functions. However, when 0 is a function of position along the duct 
perimeter, the solution for the acoustic field must be determined numerically. 
This paper is primarily concerned with analyzing the variable case. Since 
existing mathematical models cannot handle this situation, this represents 
a major contribution of this paper. 

Acoustic Equations and Boundary Conditions 
Steady state acoustic waves p>"opagating in finite ducts, such as shown 
in figure 1, are governed by the three-dimensional Helmholtz equation 

V^p + k^p = 0 (2a) 

where k = w/c is the wave number, p is the acoustic pressure, w is the 
angular frequency and c is the speed of sound. The acoustic boundary 
condition along the duct perimeter is given in the form 

^ - ik0p = 0 (2b) 

Sn 
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where 9/9n denotes the derivative in the direction of the outward normal 
from the duct wall.’ At the entrance plane, the boundary condition is 
assumed known as 

p(x,y,0) = G(x,y) (2c) 

where G(x,y) is a given input pressure. The boundary condition for the 
exit plane is expressed in the form: 


8p(x,y,d)/9z = ik0gp(x,y,d) 

(2d) 

The solution to equation (2) can be expressed in the form: 


p(x,y,z) = 1 [A e^'^n^ + B e"'''^n^] F (x,y) 
n=l ' 

(3a) 

in which 


K^= k^.x2 

n n 

(3b) 
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(3c) 
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n n nn 
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a h 

Fn(x.y) dxdy 

0 0 
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In equation (3), the characteristic function, F^, and eigenvalue, X^, satisfy 
the equation 

V^F + X^F - 0 f/|a) 

n n n ' ^ 


with the condition 


3h 



« 0 


(4b) 


along the duct perimeter. Note that as indicated in equation (3g), the 
eigenfunctions satisfying equation (4) are orthogonal. This orthogonality 
relation has been used in equation (3). 

In this paper, equation (4) will be solved by the finite element method. 
The extracted values of F^ and can then be substituted into equation (3) 
to determine the acoustic field. The numerical solution will be used to 
develop and discuss some of the attenuation characteristics of these peri- 
pherally variable liners. Equation (3) is also used to evaluate the effect 
of a finite termination impedance on the attenuation characteristics in a 
finite duct. 


Attenuation 


Before proceeding with the solution to equation (4), it is useful to 
develop an expression for the attenuation produced by the liner in terms of 
the parameters vihich have been introduced. The axial acoustic intensity at 
any axial position in the duct is 



( 5 ) 
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where is the ambient density of the medium, Re[ ] denotes the real part 
of the complex expression enclosed within the brackets and the superscript 
asterisk indicates the complex conjugate. Furthermore, the total acoustic, 
power is the integral of the acoustic intensity across the cross section 

h f. 

E = / / I-dxdy 
^ 0 0 

and the decrease in decibels of the acoustic power from z = 0 to z = Zq can 
be written as 


AdB = 10 log^oCEg/E^^) (6) 

Generally, the specific acoustic admittance, G, is chosen so as to maximize 
this attenuation. In terms of the parameters of equation (3), equation (6) 
becomes 


AdB = 10 log 


10 


^e^ It "nin \n^Q^^ 


(7a) 


in which 


and 


, i h 


''inn <"> 


ik„,z -ik^z * -ik „z 


n 


(7b) 

(7c) 


Equation (7) will allow calculation of the attenuation produced by the lining 
in the finite duct. 

Variational Formulation 

Forsythe and Wasow (ref. 9) show that equation (4) is satisfied if and 
only if the variational condition 
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fill = 0 


(8a) 


where 

II - / / I'7F- dA - ik / pF^dc (8b) 

A c 

is satisfied. In equation (0), the subscript n on F and X has been 
dropped for convenience and is the differential vector operator in the 
duct cross section. This paper will approach the solution in the duct cross 
section by using a finite element approximation to the functional II. This 
approximate functional, Tl”, is developed in the next section. 

Finite Element Formulation 

Finite elemp n^ H iscretization . - In the finite element process, the 
rectangular whose cross-section is the region (0 ^ x £ h, 0 ^ y 5 U, 
is divided into (M - 1) (fi - 1) elements as shown in Figure 2. Although 
many kinds of elements could be used in this setting, attention is restricted 
to rectangular elements. The width and height of each element are Ay and 
Ax respectively. Points at the corner of each element, called nodes, are 
designated as (I, J) where 1=1, 2, 3, . . .M and J = 1, 2, 3, . . . N. 

It is at these nodes that the values of the unknown function, F, and its 
higher derivatives will be determined. 

Shape functions. - The finite element technique expands the solution 
within each element in terms of shape functions. The shape functions at 
each node are defined in terms of local coordinates 
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Ay 

0 1 hijl £ 1 

(9b) 


The shape functions utilized in this analysis, f.|(c) and 

odd functions, respectively, which are non-zero cnly in the interval (-1, 1) 


and which have unit magnitude and slope respectively, at c = 0 as illustrated 
in figure (3). Furthermore, both f-j and f^ have zero magnitude and slope rt 
C = + 1. The functions are 



3r.^ + 2lr,)r,^ 
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Ul 1 1 
U! > 1 


(10a) 


f2(d 




(10b) 


These functions give the influence of unit magnitude and slope at the node 
in a one dimensional problem. Thus, the total approximating function may 
be formed by superimposing functions of this type. 

The approximating function for F to be used here is 

F:C(Fj 3) (11a) 
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This expansion insures that the functions F and 3F will be continuous 

t(n 

everywhere within the duct and also along its boundaries. 

The admittance along the duct pe.rimeter will be similarly expressed in 
the form 

f3 ^ 0(3^3) . (12) 

The values of the nodal coefficients in equation (12) are set to zero if 
node (I, J) does not lie along the duct wall where 3 is prescribed. Note 
that the values of the nodal coefficients in the expansion for 3 (equation 
(12)) are known whereas those in the function F, equation (11a), will be 
determined such that the functional, TT, is stationary. 

Finite element matrix equations . - The approximating function to F 
(eq. (11a)) and polynomial representation for 3 (eq. (12)) lead to an 
approximate functional TT which is a quadratic form in the unknown nodal 
coefficients Fj^ 

TT = {F}^[[S] - ik [D] tM]]{F) (13) 

where the superscr.pt T indicates the vector transpose. It is convenient 
to be consistent with the nomenclature used in structural analysis and 
refer to [S] , [D] and [M] as the acoustic stiffness, damping, and mass 
matrices respectively. [S] and [M] are square symmetric matrices whose 
coefficients are real. [D] is a complex square symmetric matrix containing 
the coefficients in equation (12), whereas (F) is an unknown vector 
containing the coefficients in the expansion given by equation (11a). [S] , 
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[M] and [D] are each of order 4MNX4MN, although many of the entries are 
zero. The unknown coefficients in {F}, (4MNX1), are ordered so that 

fpF = (f’' F^' p' 2 p22 ,11 21 

in tfll-fll, F,,. F,,, F,2_ F,2- . . F„j|} 

Requiring that the first variation of the functional vanish leads to 
a set of matrix equations of the form 


( 14 ) 


[[S] - ik[D]] {F}= A^[M] {F} 


(15) 


Equation (15) may be solved to obtain the eigenvalue X and eigenvector 


{F}. 


Finite Element Representation of Acoustic Power 
In general, the solution to equation (15) will yield many eigenvalues, 

\ and a different eigenfunction F for each eigenvalue. It is, therefore, 
convenient to refer to {F^} as the eigenvector corresponding to the eigen- 
value Now consider the integral given by equation (3g). It may be 
expressed in matrix form as: 

= {F'^MMHF^} (16a) 

likewije j 

= {f")[M]{f")* ' ( 16 b) 

Further, in order to evaluate the constant 0^ (3f)), expand the boundary 

condition G(x,y) in the form: 


G(x,y) C(g53) 

is then expressed in the matrix form as: 

= {gF[H]{F„} 




(17a) 


(17b) 
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in which 


T 1 1 

{G}' - {G 

1 1 


2 ) 12 22 11 Zi 22 

G 5 G , G , G > G ••Gw..} 
11 11 11 12 12 MN 


(17c) 


Equations (16) and (17) allow the determination of the constants and 

in the solution for the acoustic pressure (eq. (3a)) and the evaluation of 
the attenuation given by equation (7a). These expressions can then be 
utilized to evaluate the attenuation of peripherally variable liners. 


RESULTS AND DISCUSSION 

In this section, the finite element method developed in this paper is 
first compared to exact analytical solutions for uniform liners. The techni- 
que is then exercised to obtain the attenuation characteristics of infinitely 
long peripherally variable liners. Finally equation is used to evaluate 
the effects of finite duct termination on the attenuation characteristics 
of a peripherally variable liner. 

For all of these studies, certain parameters were fixed. Finite element 
re: j Its were restricted to four rows and four columns of elements (M = N = 5, 

see fig. 2) with Jl, h, and d all equal. Modes were ordered according to 

the imaginary part of the complex propagation constant. K^. The mode for 
which the imaginary part of is smallest is referred to as mode one, whereas 
the mode with the second smallest value is referred to as mode two, etc. Also, 
the bottom and two side walls of the duct are kept rigid (Bt = 33 = 62 = 0 ) and 

the eigenfunctions from the numerical and exact analysis are both normalized 

to unity at the origin. 
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Uniform Liners 


Confidence in the finite element solution may be obtained by showing it 
compares favorably with exact analytical solutions for uniform liners. 

Exact analysis . - When 3 is constant along each of the four duct 
walls (fig. 1) the solution to equation 4 is: 

F,„^(x,y) = [Cos Y^x - ^^:^S1n y^x][Cos - ^^Sin (18a) 

where the parameters satisfy the trancendental equation 


[ + k^3i3z] Sin a^z + ik(3i+ 3z) a^Cosa^il = 0 (18b) 

and the parameters Y.^j^ satisfy the same trancendental equation witha.^,3i ,3z, 

and Z replaced by Yf,»33. 3ip, and h , respectively. The eigenvalues 
are related to and Y^ by the equation 


2 2 2 


(18c) 


Watson and Lansing (ref. 10) solved the trancendental equation which has the 
general form of equation (18b) using a Newton-Raphson iterative scheme. This 
program was used in the present work to extract y^ and for the results 
involving uniform liners. These values were also used as a basis for checking 
the accuracy of the eigenvalues extracted by the finite element method. 

Comparison of the numerical and exact analysis . - The first example 
considered involved calculating eigenvalues for a hard wall duct (3 = 0). The 
lowest ten eigenvalues from equation S are listed under the heading "Exact 
Eigenvalues" in Table 1. Numerical results are also shown in the table. Good 
agreement between the finite element and exact eigenvalues is observed. The 
hard-wall modal integrals (eq. 16a) computed by the finite element method 
are also compared to exact analytical results in table 2. Good agreement is 
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again observed. 

As a second example, consider a soft wall duct in which 34 = 2.0 - 1.65i 
and the wave number is unity (k = 1). Eigenvalues extracted by the finite 
element method are compared to those extraxted by a Newton-Raphson iterative 
scheme in table 3. Good agreement is again obtained. The modal integrals 

for the soft wall duct are also compared to exact analytical values in table 4. 
Note that good comparison between the exact and analytical values is again 
obtained. 

These two examples of uniform liners illustrate the ability of the finite 
element approximations presented in this paper to give accurate and reliable 
results. Although the values are not presented, the modal integrals f n) 

were also computed for both the hard and soft wall ducts. Both the real and 

_4 

iiiuginary parts of these integrals were of the order 10 or less for each of 
the ten modes considered, indicating that the modes are orthogonal in accordance 
with equation (3g). It should also be noted that the effects of increasing 
wave number, k, have not been throughly investigated. However, it has been 
observed that at least for the first two modes, the error in both the eigen- 
values and modal integrals, / n), are consistent with those in the tables 

for 0£ k^ 10. 

Peripherally Variable Liners 

Having established confidence in the finite element solution, this 
technique is now exercised to obtain the attenuation characteristics of 
infinitely long peripherally variable liners. The effects of finite duct 
termination will also be discussed. 
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Infinitely long peripherally segmented lin&rs . - Extending the approach 
of Cremer (ref.ll), the attenuation characteristics of infinitely long 
peripheral liners are determined by investigating the attenuation rate of the 
least attenuated duct mode. The least attenuated mode in the context here 
shall mean the mode with the smallest attenuation rate. Cremer (ref.ll) has 
determined thr 4 ti for a two-dimensional infinitely long duct with a hard 
surface on the lower wall, the attenuation rate of the least attenuated duct 
mode can be maximized by choosing the admittance of the upper wall to be: 

3(,(kh) = (2.06 - 1.65i)/kh (19) 

A liner with a wall admittance of 3_(kh) is referred to as a Cremer liner and 
is said to be tuned at Kq if the admittance of the upper wall is given by 
g^(kgh). This relation has been utilized to determine admittance values for 
the peripherally segmented liners. 

In the examples to follow, the bottom and two side walls of the duct in 
figure 1 have been kept rigid (3i = 3a = 3z =0). Thus the admittance 
(see fig. 1) of the upper wall is allowed to vary in a stepwise or continuous 
fashion. This is at variance with the mathematical models which have been 
constructed to date which a'p restricted to axial segmentation of the duct 
liner (segmentation along the Z-axis). 

In order to simplify discussions of peripherally segmented liners, a 
segmented liner will be referred to as a k^-kg-k^-k^ liner. This notation 
implies that: 
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3c(l<2h ) 

3c(k3h ) 


0 £ y < jl' 

•n 

•< 2 

1 < y £ j. 

2 ■• 

i.£y < 1 

4 


( 20 ) 


where 3-(kh) is the Cremer admittance given by equation 19. Also, to 
avoid referring to a Cremer segment tuned at infinity, the notation 
- 0 - kg - k^ will denote that the portion of 3 ^ corresponding 
to the tuning frequency kg is a hard wall. 

Numerical studies were conducted to evaluate the effects of 
shorteningthe width of an initially full width Cremer liner. Figure 4 
illustrates this effect. Note that once a small portion of the liner is 
taken away, the attenuation rate near the tuning frequency decreases 
substantially. On first thought, one might expect that by taking away 
one quarter of the liner, the attenuation rate near the tuning frequency 
would be reduced by about one quarter. However, in fact the attenuation 
rate is reduced by about one half. 
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In figure 5, it can be seen that the attenuation rate of a peripherally 
segmented liner may also be changed simply by changing the positions of 
the liner segments. Note, that although the attenuation rate near the 
tuning frequency of the liner remains the same, the attenuation curve in 
the mid-frequency range is much broader for a 0-2-0-2 or 0-2-2-0 liner than 
for a 2-2-0-0 liner. 

In figure 6, the attenuation characteristics of peripherally segmented 
liners having Cremer segments tuned at several frequencies are shown. The 
first segment of all liners is a Cremer liner tuned at h = 2. It can be 
seen that the shape of the attenuation curve can be altered, although not 
substantially, by the variation of the tuning frequencies. Essentially, 
the curve is dominated by the lowest tuning frequency utilized. 


Infinitely long continuous peripheral liners . - The finite element 
analysis presented in this paper can also be employed to investigate a periphe 
ral liner with continuous variation. 

Figure 7 shows the attenuation curves of four peripheral liners with 
continuous linear variation of the admittance. The admittance Pi, of each 
liner is given by: 


where 


3/1 = An + A, y 


- 3^(2. 0 h) 


^1 " ^C'*^max 


(21a) 


(21b) 

(21c) 
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The value of Is different for each liner. Note that for each liner 

the maximum attenuation rate occurs at h = 2. However, the value of this 
maximum attenuation decreases with increasing k,„^h. There is also some 

IIiua 

broadening of the attenuation curve at the midfrequencies. 

In figure 8, the variation of the real part of the admittance of three 

different types of liners are compared. The segmented liner is a 2-4-6-8 

liner, while the liner with the liner variation has the form of equation 21 

with k„,„h=8. The admittance for the cubic liner was expressed in the form: 
max 

34(y) = Aq + A^y + Agy + (22) 

where the coefficients in equation 22 v/ere chosen to that, aty = 0, y = 1/3, 
y = 2/3, and y = 1 , the values of were 6^(2. 0 h), 8^(4. 0 h), 3c(6.0 h), 
and 8^(8. 0 h) respectively. The variation of the imaginary part of the 
admittance of these three liners will be similar to that shown in figure 8. 

Figure 9 shows the attenuation rates for each of the three liners. The 
attenuation rates for the cubic and four-segment liner are not significantly 
different although the attenuation curve for the cubic liner is slightly 
broader than that for the segmented liner. The liner with linear variation 
has a greater attenuation rate than the other two liners in the low to 
midfrequency range. 

Effects of a Termination Impedance on Liner Attenuation 
When the duct is of finite length, non-negligible acoustic energy may 
be carried to the termination by other modes than the least attenuated mode. 

As a result, the attenuation produced by the liner must be computed from 
equation 7a, which includes the effect of the termination as well as the 
interaction between, the different modes. At the present time, the computer 


program utilized in this study 1s unable to compute this equation exactly 

due to storage limitations. The eigenfunctions , , and modal integrals 

and must all be retained. In addition, since the subroutine that 
nn mn 

solves the eigenvalue equations (eq. 15) destroys the mass matrix, the 
coefficients in this matrix must be stored in another matrix, thus requiring 
further core. 

However, if it is assumed that the source is such that all modes have 
equal input energy and are in phase and that coupling between the modes can 
be neglected = 0, m n), then the attenuation in the finite duct is 
proportional to the sum of the attenuations of each of the individual modes. 
Equation 7a reduces in the case of a single mode to: 


AdB = 10 Log,p Rele n -c„e n + c^e n c„c^e n j (jjj 

Re[l - + c^c* ] 

Equation 23 is the expression for the attenuation of a single mode in a 
finite duct. For infinitely long ducts, c^^ = 0. However, the parameter c^ 
is not necessarily zero for each mode in the finite duct and thus will 
affect the value of the attenuation for the mode. Note that c^(eq. 3h) 
contains the effect of the termination admittance as well as the duct length, d. 


In the remainder of this section, the effects of a p^c termination 

impedance (3 =1) upon the attenuation of a single mode in a finite duct 

will be investigated. Note that since the length of the duct is unity, 

the expression on the right-hand side of equation 23 may also be referred 

to as the attenuation "rate". When c 0, the mode for which the right- 
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hand side of equation 23 is smallest will be referred to as the least 
energy mode. Hov/ever, when = 0, the mode for which the right-hand 
side of equation 23 is smallest will be referred to as the least attenuated 
mode , 

The attenuation rates for a 2-4-6-S liner and a 2-2-8-8 liner are 
given in figures 10 and 11, respectively. In each figure, the attenuation 
curves for the least energy mode, the least attenuated and the least 
attenuated mode with reflections (c^ f 0) are plotted. The asterisk on 
the attenuation curve for the least energy mode denotes that, at that parti- 
cular value of kh, the least energy mode is the same as the least attenuated 
mode. Note that figures 10 and 11 display the same characteristics. The 
attenuation rate for the least energy mode is lesi; than that of the least 
attenuated mode with reflections at the lower frequencies. However, each 
of the three attenuation curves are identical at the higher frequencies. 

In fact, the least attenuated mode and the least energy mode are the same 
mode at the higher frequencies. 

The results shown on figures 10 and 11 indicate that the termination of 
the duct must be taken into account in liner optimization studies at the 
lower frequencies. At the higher frequencies, the attenuation rates are 
identical indicating that the effects of the termination can be neglected 
at these frequencies. One should understand however, that these conclusions 
are based on the assumption that the total attenuation rate in the finite 
duct is proportional to the attenuation rate of the individual niodes. 




CONCLUDING REMARKS 


Based on the results of this work, the following conclusions are drawn: 

1. The attenuation rate for a peripherally segmented liner drops 
significantly if a small portion of the liner is removed. 

2. One can alter the attenuation characteristics of a peripherally 
segmented liner by placing the liner segments in different 
positions along the periphery. 

3. Based upon the attenuation curves considered here, it appears 
that the primary effect of peripheral variation is a broadening 
of the attenuation characteristics for the liner. 

4. Effects of the duct termination must be taken into account at 
the lower frequencies, when optimizing peripherally variable 
liners. 

Conclusions 1, 2, and 3 were based upon infinite ducts, whereas conclusion 4 
employed the assumption of equal partition of energy and neglecting coupling 
for the modes. Although these restrictions are somewhat removed from reality, 
they do provide a means for investigating the attenuation characteristics of 
peripherally variable liners. 

Future Work 

There is additional work to be done with regard to the computer program 
developed in conjunction with the finite element theory. In order for the 
computer program to model more realistic ducts, such as finite length ducts 
with realistic sources, it is suggested that equation 15 be solved out of core. 
Both matrices are block tridiagonal, symmetric, and in addition [M] and [S] are 
real. Taking advantage of these special features should allow one to incor- 
porate thousands of degrees of freedom in this present analysis. It is further 
suggested that the computer program be combined with an optimization program so 
as to determine if peripherally variable liners can achieve more attenuation 
than the best uniform liners. 
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TABLE 1 


COMPARISON OF EXACT AND FINITE ELEMENT EI'iENVA.JES 
FOR A HARO WALL OUCT (6^-0) 


'TCDE no., n 

EXACT EIGENVALUE 

eigenvalue COMPUTED I 
FROM FINITE ELEMENT 1 
METHOD 

1 

0.000 

0.000 

2 

3.142 

3.142 

3 

3.142 

3.142 

4 

4.443 

4.443 

5 

6.2B3 

6.284 

6 

6.283 

6.284 

7 

7.025 

7.025 

8 

7.025 

7.027 

9 

8.886 

8.887 

10 

9.425 

1 

9.434 


TABLE 2 

COMPARISON OF EXACT MODAL INTEGRALS WITH THOSE EXTRACTED 
BY THE FINITE ELEMENT METHOD FOR A HARD WALL DUCT 


MODE NO., n 

EXACT V 

nn 

VALUE OF 'V COMPUTED 
nn 

FROM FINITE ELEMENT 
METHOD 

1 

1.000 

i 1.000 

2 

.500 

I .500 

3 

.500 

i .501 

f 

4 

. .250 

! .250 

5 

.500 

1 .500 

6 

.500 

1 .495 

7 

.250 

.246 

8 

.250 

.247 

9 

.250 

.245 

10 

.500 

.480 1 





TABLE 3 


COMPARISON OF EIGENVALUES EXTRACTED BY THE FINITE ELEMENT 
METHOD WITH THOSE OBTAINED BY THE NEWTON-RAPHSON 

ITERATIVE SCHEME FOR k - 1. 6 ^-B 3 -B 2 ■ 0 and 2.0-1.651 


MODE NO., n 

EIGENVALUE COMPUTED FROM 
NEWTON-RAPHSON ITERATIVE 
SCHEME 

EIGENVALUE COMPUTED 
FROM FINITE ELEMENT 
METHOD 

1 

1.835 - 1.2191 


2 

2.342 - .987! 


3 

3.487 - .6411 


4 

3.840 - .6021 

3.844 - .6021 

5 

6.012 - .3361 

6.013 - .3371 

6 

6.440 - .3471 

6.432 - .3381 

7 

6.542 - .3481 

6.656 - .3511 

8 

6.782 - .2981 

6.783 - .2931 

9 

8.693 - .2331 

8.696 - .2321 

10 

9.250 - .2171 

9.256 - .2181 


TABLE 4 

COMPARISON OF THE EXACT MODAL INTEGRALS WITH THOSE COMPUTED BY 
THE FINITE ELEMENT METHOD FOR THE DUCT 
OF TABLE 3 
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Figure 1. - Schematic of three-dimensional finite duct with peripheraliy variable liner. 
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Figure 2. - FIniie elenient c* I sc ret ion of (^uct cross-section. 
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Figure 6. - Attenuation rate of a segmented liner tuned at several frequencies. 
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FOUR SEGMENT LINER 
LINER WITH LINEAR VARIATION 



Figurs 5. * Aoirniiancy vai ialion far a four segraent, linear ar 
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LEAST AHENUATEO MODE WITH REFLECTIONS 

LEAST ENERGY MODE 

MODE CORRESPOND TO LEAST AHENUATED MODE 



LEAST ATTENUATED MODE 

LEAST ATTENUATED MODE WITH REFLECTIONS 

LEAST ENERGY MODE 
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Figure 11.- Attenuation rate in a finite and infinite duct 
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10 AbtDOC' 

Using multimodal analysis, a variational finite element method 
is presented for analyzing sound attenuation in a three-dimensional 
finite duct with a peripherally variable liner in the absence of flow. 

A rectangular element, with cubic shape functions, is employed. 

Excellent comparison between exact results and the finite element motho( 
is obtained in cases where an analytical solution exists. This study 
indicates that, once a small portion of a peripheral liner is removed, 
the attenuation rate near the frequency where maximum attenuation 
occurs drops significantly. Also, it was observed that the positioning 
of the liner segments affects the attenuation characteristics of the 
liner and that effects of the duct termination are important in the 
low frequency ranges. In general, the results indicate that the main ^ 
effect of peripheral variation of the liner is a broadening of the 
attenuation characteristics in the mid-frequency range. Finally, it is; 
concluded that, due to matrix size limitations of the presently avail- ' 
able computer program, the eigenvalue equations should bo solved out 
of core in order to handle realistic sources. 
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